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ABSTRACT 

Let  X  be  a  given  a  set  of  n  circular  and  straight  segments  in 
the  plane  where  two  segments  may  only  intersect  only  at  their 
endpoints.  We  introduce  a  new  technique  that  computes  the 
Voronoi  diagram  of  X  in  O(nlogn)  time.  This  result 
simultaneously  improves  on  several  previous  algorithms  for 
special  cases  of  the  problem.  The  new  algorithm  is  relatively 
simple  compared  to  previous  algorithms,  a  fact  important  for 
the  numerous  practical  applications  of  the  Voronoi  diagram. 
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1.   Introduction 

The  ubiquitous  Voronoi  diagram  has  been  extensively  studied  in  areas  such  as  biology, 
solid-state  physics,  pattern  recognition,  geography,  stock-cutting,  wire  layout,  geometric 
optimization,  facilities  location,  computer  graphics,  and  robotics  (see  [Kir79],[Dry79],  or 
[Lee82]  for  extensive  references).  In  some  literature,  the  alternative  terminology  of  Thiessen 
or  DirichUt  tesselation  is  more  common.  In  this  paper  we  are  concerned  with  the  Voronoi 
diagram  for  a  set  of  planar  objects  under  the  Euclidean  metric.  When  restricted  to  the 
interior  of  a  simple  polygon,  this  diagram  is  known  as  the  medial  axis  or  internal  skeleton  of 
the  polygon.  Many  variations  of  the  problem  studied  here  have  been  investigated.  To 
illustrate  the  range  of  possibilities,  we  mention  the  following: 

(a)  Using  the  general  L  metric  instead  of  the  usual  Euclidean  metric.  An  unusual  metric 
that  arises  in  computational  fluid  dynamics  [WPB]  is  the  problem  of  computing  the 
Voronoi  diagram  for  a  set  of  points  under  the  following  metric: 

D(p,  q)  =  Tnin{d(p,  q+r)  :  r  ^  Z'}, 

where  Z  are  the  integers,  and  d(p,  q)  the  Euclidean  distance.  We  may  interprete  this  as 
computing  the  Voronoi  diagram  on  the  torus.'  This  problem  can  be  linearly  reduced  to 
the  standard  Voronoi  diagram  for  a  set  of  points. 

(b)  Assigning  weights  to  points  [AES3]. 

(c)  Generalization  to  higher  dimensional  spaces.  An  unusual  space  is  R^  x  S'  where  R  is  the 
real  line  and  S'  the  unit  drcle:  the  Voronoi  diagram  here  is  used  for  planning  the  motion 
of  a  line  segment  ([OSY83],[OSY84a],[OSY84b]). 


'  The  points  on  the  torus  rqjresent  (moving)  marlcen  for  sdving  Navier-Stokes  equations. 


Our  interest  in  the  problem  arose  from  robctia:  m  [OY82]  we  show  that  planning  the 
motion  of  a  disc  amidst  polygonal  objects  can  be  reduced  to  searching  in  the  Voronoi  diagram. 
This  result  dearly  extends  to  the  case  where  the  obstacles  are  bounded  by  line  segments  or 
circular  arcs.  (As  an  example  of  the  use  of  circular  arcs,  the  mobile  robot  in  [MorSl] 
approximates  itself  and  the  obstacles  by  discs;  an  0(n*)  time  algorithm  was  implemented 
there.)  The  importance  of  the  case  of  a  disc  is  the  algorithm's  extreme  efficiency  compared  to 
other  even  slightly  more  complicated  shapes  (eg.  [OSY83]).  Indeed,  computing  the  Voronoi 
diagram  can  be  regarded  as  a  preprocessing  cost  in  which  case  the  actual  path  planning 
becomes  a  linear  time  process.  Thus  an  algorithm  for  moving  a  disc  can  serve  as  an  important 
initial  heuristic  in  general  motion-planning  algorithms.  But  until  the  availability  of  an  easily 
implementable  0(n\opi)  algorithm  for  computing  the  Voronoi  diagram,  this  importance  may 
remain  mostly  theoretical.  The  results  of  the  present  paper  will  hopefully  remedy  this 
situation.  (See  also  [Siia83]  for  another  robotics  application  of  Voronoi  diagrams:  here  we 
want  to  detect  if  any  two  differently  colored  circles  intersect  from  a  set  of  colored  circles.) 

Previous  Work.  Before  the  advent  of  computational  geometry,  a  number  of  algorithms  for 
various  cases  of  the  problem  considered  in  this  pjaper  were  proposed  (mostly  running  in  time 
fl(n'^).  Here  we  wiU  review  those  results  that  are  asymptotically  efficient  and  which  rely  on 
the  techniques  of  computational  geometry. 

(a)  The  first  algorithm  in  this  genre  is  due  to  Shamos  and  Hocy  [SH75]  who  gave  an 
O(nlogn)  algorithm  for  the  Voronoi  diagram  for  a  set  of  points.  We  might  add  that 
algorithms  for  the  Voronoi  diagram  of  point  sets  in  higher  dimensions  follows  indirectly 
from  the  work  of  Seidel  on  the  convex  hull  [SeiSl]  , 

(b)  In  the  thesis  of  Drysdale  [Dry 79],  the  problem  of  the  Voronoi  diagram  for  a  set  of 
disjoint  polygonal  and  circular  objects  was  studied.  He  described  and  implemented  an 
0{nc  '°8'')  algorithm.  This  bound  is  subquadratic  but  since  it  is  n(nlog^n)  for  any  k, 
Drysdale  also  posed  the  problem  (solved  in  this  paper)  of  an  0{n\ogn)  solution. 

(c)  Subsequently,  Dr>'sdale  and  Lee  improved  the  bound  in  (b)  to  0(n!og-n)  [LD81]. 

(d)  At  the  1979  FOCS,  Kirkpatrick  [Kir79]  outlined  an  C>(/ilogn)  solution.  But  the 
technique  is  complicated  enough  that  the  details  still  remain  somewhat  unsettled 
([Kir84],[Sha84]).  Our  techniques  is  simpler  than  Kirkpatrick's  although  his  ideas 
('spokes'  and  the  use  of  minimum  spanning  tree)  have  independent  interest. 

(e)  When  restricted  to  the  problem  of  computing  the  medial  axis  of  a  simple  polygon, 
[Lee82]  presented  an  O(nlogn)  algorithm  (improving  an  earlier  one  in  [P*re77]). 

(f)  aiarir  [Sha83]  describes  an  O(nlog-n)  algorithm  for  circles  that  may  intersect.  Note  that 
this  improves  (c)  for  the  case  of  drdes  since  the  solution  in  (c)  assumes  that  the  drdes 
are  disjoint.  We  remark  that  the  connection  between  Sharir's  problem  to  the  setting 
assumed  in  this  paper  is  not  direct.  At  the  end  of  the  paper  we  will  show  how  a  variant 
of  our  technique  improves  his  bound  to  O(nlogn). 


We  refer  to  a  recent  review  [SL85]  of  most  of  the  preceding  results  as  well  as  the 
techniques  and  generalizations  known. 

Technical  discussion.  All  the  above  algorithms  use  the  divide-and-conquer  paradigm:  To 
compute  the  Voronoi  diagram  of  X,  divide  X  into  equal  subsets  X^  and  X^,  recursively 
compute  their  Voronoi  diagrams  and  then  'merge'  the  result.  The  merging  is  essentially 
defined  by  a  certain  'merge  curve'  C.  To  obtain  an  O(nlogn)  algorithm,  the  goal  is  to 
compute  C  in  linear  time.  In  the  Shamos-Hoey  algorithm,  the  merge  curve  is  a  connected  set 
(we  can  view  this  as  a  kind  of  'separability'  property  of  the  two  sets  X^  and  X/().  The  work  of 
Drysdale  and  Lee  attempts  to  recover  this  separability  property.  As  they  reported,  finding 
such  a  separability  property  that  is  computationally  simple  remained  elusive  despite 
considerable  effort.  Accepting  the  fact  that  C  will  have  many  connected  components  in 
general,  the  technical  issue  is  to  t"md  at  least  one  pomt  (called  a  'starter')  in  eacli  compwnent 
of  C.  The  innovation  of  Kirkpatrick  is  to  show  that  no  notion  of  separability  is  needed  (ie.  X^ 
and  X^  can  be  arbitrary).  His  idea  is  to  subdivide  each  Voronoi  cell  (by  introducing  'spokes') 
into  simpler  subcells,  and  to  use  the  fact  that  a  certain  minimum  spanning  tree  of  X  intersects 
the  Voronoi  edges  and  spokes  of  X^  and  X^  in  a  fashion  that  allows  one  to  find  the  starters. 
This  idea  appeared  again  in  Sharir's  work  on  intersection  circles.  In  some  sense,  our  new  idea 
is  to  reintroduce  the  separability  condition  in  a  radical  way  ('by  simply  imposing  it'). 

Underlying  technique  of  this  paper.     We  now  give  the  basic  idea  of  this  paper.    Recall  from 

the  abstract  that  the  problem  we  solve  is  the  following: 

Given  a  set  X  consisting  of  n  straight  or  circular  segments  (possibly  degenerated  to 
points),  where  the  segments  do  not  intersect  except  at  their  endpoints,  compute  their 
Voronoi  diagram  Vor{X). 

This  formulation  subsumes  (generalizes)  all  the  work  expliddy  mentioned  above  since 
polygons  and  circles  can  be  decomposed  into  circular  and  straight  segments.  It  will  be  shown 
in  the  next  section  that  the  Voronoi  diagram  is  composed  of  straight,  parabolic,  hyperbolic  or 
elliptic  curves.  Since  all  and  only  conies  can  appear  in  the  diagram,  our  problem  is  a  very 
natural  level  of  generalization  of  the  original  problem  for  points. 

If  there  are  m  distinct  endpoints  among  the  segments  of  X,  we  introduce  m  +  1  vertical 
lines  such  that  every  endpoint  of  a  segment  of  X  lies  between  a  unique  pair  of  adjacent 
vertical  lines.  The  region  between  any  pair  of  vertical  lines  is  called  a  slab.  In  stage  0,  for 
each  slab  between  a  pair  of  adjacent  vertical  lines,  we  compute  the  Voronoi  diagram  of  the 
restriction  of  the  segments  of  X  to  the  slab.  In  stage  1,  we  'merge'  the  Voronoi  diagrams  of 
adjacent  slabs  computed  in  stage  0.  In  general,  at  stage  i  +  l  we  merge  pairs  of  adjacent  slabs 
from  stage  i.  In  log  m  stages,  we  would  essentially  have  computed  the  the  Voronoi  diagram 
of  X.  Note  that  in  stage  i  we  compute  the  Voronoi  diagram  of  slabs  that  contains  2' 
endpoints  of  X.  The  obvious  implementation  of  this  algorithm  takes  CL(n~)  time,  simply 
because  in  the  first  stage,  merging  each  pair  of  slabs  can  take  up  to  linear  time. 


We  overcome  this  problon  by  computing  only  the  'essential'  part  of  the  Voronoi 
diagram  of  a  slab,  where  this  essential  part  has  size  only  0{k)  if  the  slab  contains  k  cndpoints 
of  X.  Furthermore,  merging  two  slabs  that  collectively  contain  it  endpoints  takes  only  0{k) 
work.  (Note:  actually,  k  should  be  the  number  of  segments  with  endpoints  in  the  slab.)  This 
implies  that  each  stage  takes  linear  time  and  our  stated  time  bound  follows. 

In  the  rest  of  this  paper,  we  first  present  the  basic  definitions  in  section  2.  Section  3 
describes  some  simple  properties  of  the  process  of  moving  along  the  Voronoi  diagram  of 
objects.  Section  4  gives  the  merging  process  which  is  the  heart  of  the  algorithm.  Section  5 
and  6,  respectively,  prove  the  termination  and  the  correctness  of  the  merge  process.-  After  the 
overall  algorithm  is  given  in  section  7,  we  analyze  its  complexity  in  section  8.  In  section  9,  we 
show  the  techniques  in  this  paper  lead  to  an  0(/ilogn)  solution  to  Siarir's  variant  of  the 
Voronoi  diagram  for  n  intersecting  circles.   We  conclude  in  section  10. 

2.   Preliminaries 

Following  [Kir79]  we  take  our  primitive  objects  to  be  points,  opjen  line  segments  and 
open  circular  arcs.  It  is  important  to  remark  that  this  'expedient'  of  Kirkpatrick  is  actually  a 
crucial  insight  that  solves  several  tedinicaJ  problems  faced  when  trying  to  generalize  the 
original  definition  of  Voronoi  diagrams  for  points:  sec  [Dry 79]  for  an  extensive  discussion  of 
these  issues.  For  simplicity,  we  restrict  the  line  segments  to  be  finite  and  the  circular  arcs  to 
be  less  than  a  semi-drde.  Typically  the  line  segments  come  from  a  set  of  polygonal  objects 
and  the  arcs  come  from  cirdes.  Motivated  by  robotics  applications  (eg.  [OY82]),  the  open 
line  segments  and  arcs  will  be  called  walls  and  the  points  will  be  called  conurs.  An  object  is 
either  a  wall  or  a  comer.  (Note:  a  comer  is  treated  as  a  point  or  a  singleton  set  as 
convenient.)  A  set  X  of  objects  is  said  to  be  proper  if  (a)  they  are  pairwise  disjoint  and  (b)  for 
each  wall  in  X  its  endpoints  are  comers  in  X.  Note  that  we  allow  isolated  comers  in  a  proper 
set.  Hereafter,  let  X  denote  a  proper  set  of  objects.  We  remark  that  the  assumption  (b)  in 
the  definition  of  'propemess'  is  a  technical  convenience  which  ensures  that  each  Voronoi  edge 
is  part  of  a  unique  conic  (see  below). 

Some  of  the  following  definitions  are  fairly  standard:  The  projection  of  a  point  p  onto  an 
object  s  is  the  point  <?  in  the  dosure  7  of  j  such  that  the  Euclidean  distance  d{p,  q)  is 
minimized  [LD81].  Note  that  this  definition  is  well-defined  except  when  j  is  a  circular  wall 
and  p  is  at  the  center  of  the  drde  containing  s.  The  distance  d(p,  s)  between  p  and  s  is 
defined  as  inf{d{p,  q)  :  q  i  s).  Define  the  clearance  of  an  arbitrary  point  p  with  respect  to  X 
to  be  the  minimum  of  d{p,  s)  where  s  i  X.  Denote  this  by  Clearance -i^)  or  simply 
Cleararxceip).  In  our  proofs,  it  is  often  convenient  to  refer  to  the  drde  centered  at  p  with 
radius  Clearancexip)'-  we  call  it  the  clearance  circU  at  p  (with  respect  to  X). 


'  As  Sharif  [Sha83]  panted  out,  the  con-ectness  and  termination  of  the  Lee-Drysdile  scanning  procedure 
(wtiich  our  algorithm  depends  on)  is  not  at  all  ob\'ious. 


We  consider  two  ways  to  define  the  Voronoi  diagram.  K  we  regard  X  as  a  set  of  points 
obtained  as  the  union  of  the  objects  in  X  then  we  have  a  very  simple  'intrinsic'  definition  as  in 
[OY82]  or  [OSY84a].  Precisely,  if  \JX  is  the  set  of  points  in  objects  of  X,  then  the  inn-insic 
Voronoi  diagram  of  X  is  the  set  of  points  p  such  that  the  intersection  of  the  clearance  circle  at 
p  with  [JX  is  a  disconnected  set.'  The  defmition  of  the  intrinsic  diagram  serves  to  motivate 
the  next  definition  of  the  Voronoi  diagram.  It  is  essentially  due  to  [Kir 79].  Say  a  point  p  is 
cbse  to  an  object  jinXifforaIlt>0  there  is  a  point  (?  in  the  t-neighborhood  of  p  such  that 
(i)  Clearance{q)  =  d{q,  s)  and  (ii)  the  projection  of  q  onto  j  is  actually  in  s  (rather  than  in 
J-s).  It  is  not  hard  to  observe  that  p  is  dose  to  s  iff  CUarance(p)  =  d(p,  s)  and  either  (a) 
the  projection  of  p  onto  s  is  indeed  in  j,  or  (b)  if  j  is  a  wall  and  p  is  on  the  normal  through 
an  endpoint  of  s . 

Definition.  The  Voronoi  diagram  Vor{X)  of  X  is  the  set  of  points  p  such  that  there  exist 
two  objects  s,  and  Ji  that  are  dose  p. 

By  the  last  observation,  Vor(X)  is  simply  the  intrinsic  Voronoi  diagram  of  X  augmented  by 
additional  line  segments  lying  along  the  normals  to  the  endpoints  of  walls  (case  (b)  above). 
For  computational  purposes  and  for  most  of  this  paper,  we  prefer  to  use  Vor{X)\  but  we  will 
have  a  few  occasions  to  consider  the  intrinsic  diagram  of  X.  These  are  the  cases  where  the 
objects  in  X  may  arbitrarily  intersect,  and  the  resulting  Vor(X)  becomes  problematic  as 
discussed  by  Drysdale. 

For  any  pair  of  disjoint  objects  s,  s' ,  the  (s,  s')-bisector  is  the  Voronoi  diagram  of  the 
set  (j,  s'}.  In  this  paper,  we  only  consider  (j,  j')-bisectors  where  s  and  s'  are  objects  from  X. 
The  pTOpemess  of  X  implies  that  such  a  bisector  is  a  simple  curve  that  divides  the  plane  into 
two  Lnfmite  regions.  But  in  general,  the  bisector  may  contain  branch  points  as  illustrated  by 
the  following:  If  s  and  s'  are  two  straight  walls  such  that  one  endpoint  <7  of  j  is  in  the  interior 
of  s'  then  the  bisector  consists  of  three  branches  emanating  from  q.  (This  situation  is 
exduded  by  the  propemess  of  X  since  q  would  have  to  be  an  object  in  X  but  <?  and  s'  are  not 
disjoint.) 

Let  us  briefly  note  the  types  of  bisectors:  in  case  the  objects  are  comers  or  straight  walls, 
the  bisector  is  familiar  from  previous  work.  For  instance,  ii  s,  s'  are  both  comers  then  the 
bisector  is  a  line;  if  j'  is  a  straight  wall  and  s  is  its  endpoint  then  the  bisector  is  the  line 
through  s  and  normal  to  s' ;  if  they  are  both  straight  walls  then  the  bisector  is  a  curve  that 
may  be  composed  of  up  to  7  sections  of  straight  or  parabolic  lines  (note:  this  case  is  essentially 
figure  2  if  we  ignore  the  endpoints  of  the  two  walls  and  the  related  Voronoi  edges).  Next  wc 
illustrate  the  basic  types  of  interactions  involving  circular  arcs.  It  is  enough  to  use  infinite 
lines  and  full  cirdes  instead  of  line  segments  and  arcs.  We  also  allow  these  cirdes  and  lines  to 
freely  intersect  and  hence  we  will  use  the  intrinsic  Voronoi  diagram  here.    It  is  then  easy  to 


'  This  elegant  definition  of  VoronQi  diagraits  is  an  alternative  way  to  overcome  the  atxA'e-mentioned 
difficulties  in  defining  Voronoi  diagrams. 


verify  the  following: 

(1)  Two  non-intersecring  circles:  the  bisector  is  a  hyperbola  (unless  the  two  radii  are  equal, 
in  which  case  we  have  a  straight  line). 

(2)  Two  intersecting  circles,  neither  containing  the  other,  the  bisector  is  the  union  of  an  ellipse 
and  a  hyperbola,  both  passing  through  the  two  intersecting  points  of  the  circles.  See 
figure. 

(3)  Two  circles,  one  contained  in  the  other:  the  bisector  is  an  ellipse  with  fed  the  two  centers 
of  the  circles.  The  ellipse  separates  the  two  circles.  In  the  degenerate  case  where  the 
two  circles  touch,  there  is  an  additional  ray  from  the  center  of  the  contained  circle 
directed  away  from  the  other  center. 

(4)  A  circle  of  radius  r  and  a  non-intersecting  line  L:  a  parabola  with  focus  the  center  of  the 
drcle  and  directrix  a  line  parallel  to  L  and  at  distance  r  from  L . 

(5)  A  circle  of  radius  r  and  an  intersecting  line  L:  two  parabolas  both  passing  through  tlie 
two  intersection  points  of  L  and  circle.  The  directrices  of  the  two  parabolas  are  the  two 
lines  parallel  to  L  and  at  distance  r  from  L . 

(6)  A  circle  and  a  point  outside  the  circle:  a  hyperbola. 

(7)  A  circle  and  a  point  inside  the  circle:  an  ellipse  inside  the  drde.  (Note:  in  the  analysis, 
we  can  essentially  treat  cases  (6)  and  (7)  as  degeneracies  of  (1)  and  (3),  respectively.) 


Figure  1.  Bisector  of  two  intersecting  circles 

Remark:  The  reader  familiar  with  the  work  of  [Sha83]  will  note  that  our  definition  of 
Voronoi  diagrams  (when  restricted  to  full  circles)  differs  from  Sharir's  (which  is  similar  to  that 
in  [LX)81]).  Sharir  works  with  drcles  and  defines  the  distance  D(p,  C)  from  a  point  p  to  a 
drcle  C  of  radius  r  centered  at  <7  as  d{p,q)-r.  So  distance  could  be  negative  and  the 
diagram  for  a  set  of  cirdcs  defined  by  Sharir  has  no  elliptic  curves.    But  it  is  easy  to  see  that 


by  removing  all  the  elliptic  in  our  diagram,  we  obtain  Sharir's.   The  disadvantage  of  negative 
distances  is  that  it  is  not  easy  to  generalize  to  droilar  arcs. 

The  Voronoi  diagram  can  be  decomposed  into  Vororwi  edges  where  each  edge  e  is 
determined  by  a  pair  of  objects  j,  j' ;  in  this  case  e  is  called  an  (s,  s')-edge.  It  can  be  shown, 
by  adapting  standard  proofs,  that  the  number  of  Voronoi  edges  is  linear.  It  is  clear  that  the 
(j,  j')-edge  in  Vor{X)  is  a  portion  of  the  (j,  i')-bisector.  Because  X  is  proper,  each  edge  is  a 
segment  of  a  unique  conic  (rather  than  a  union  of  such  segments).  The  endpoints  of  Voronoi 
edges  are  called  Voronoi  vertices.  The  set  of  points  in  the  plane  that  are  neither  on  the  Vor{X) 
nor  in  any  object  of  X  is  partitioned  into  connected  components  called  Voronoi  cells.  For 
example,  there  are  8  cells  in  the  next  figure. 


f,*- 


-4r. 


Figure  2.  Voronoi  cells 


Each  cell  is  associated  with  an  object  s  where  for  all  points  p  in  the  cell,  p  is  dose  to  j  in  the 
above  sense.  Conversely,  a  comer  (resp.  wall)  is  associated  with  at  most  one  (resp.  exactly 
two)  cells.   A  cell  associated  with  an  object  s  is  called  a  s-cell. 

The  remaining  definitions  have  been  invented  for  the  technique  in  this  paper.  Let 
m  5  n  be  the  number  of  comers  in  X.  As  in  the  introduction,  let  us  introduce  m  +  \  vertical 
lines  called  separators  such  that  each  comer  is  between  a  unique  pair  of  adjacent  separators. 
The  strip  between  any  two  separators  is  called  a  slab.  We  may  assume  that  a  circular  wall  has 
the  property  that  any  vertical  line  intersects  it  at  most  once. 


Let  5  be  a  slab.  A  wall  is  long  with  respect  to  5  if  it  intersects  both  of  the  separators 
that  bound  S.  The  set  of  long  walls  partitions  the  slab  into  closed  regions  that  we  call  quads, 
so-called  because  these  have  four  sides  when  they  are  finite  regions.  Thus  if  a  slab  has  no 
long  walls  then  the  whole  slab  is  the  quad;  otherwise,  all  but  two  of  the  quads  are  finite.  A 
quad  is  said  to  be  active  if  there  is  at  least  one  comer  in  it;  otherwise  it  is  inactive.  Let  I  be  a 
separator.  A  crossing  (or  s -crossing)  is  the  intersection  of  a  separator  with  a  wall  s.  If  Z,  has 
k  ^  0  aossings,  then  L  is  divided  m  k  +  1  segments  called  windows.  A  window  is  active  with 
respect  to  5  if  it  is  part  of  the  boundary  of  an  active  quad  of  S. 
Example: 


Figure  3.  A  slab  between  two  separators  (the  dotted  lines) 


In  the  figure,  we  have  a  slab  with  4  (not  5!)  quads  of  which  the  top  two  quads  are  active.   The 
second  quad  from  the  top  has  3  windows  on  the  left  boundary  and  2  on  the  right. 

Let  2  be  a  quad.  A  Q -object  j  is  an  object  that  is  either  a  comer  corresponding  to  a 
crossing  at  a  vertical  boundary  of  2  or  an  object  of  the  form  s  -  s'  nint(Q)  where  s'  is  an 
object  of  X  and  int(Q)  is  the  interior  of  Q.  We  also  refer  to  a  Q-object  as  a  Q-wall  or  Q- 
comer  as  the  case  may  be.  For  instance,  if  j'  is  a  long  wall  that  defines  the  upper  or  lower 
boundary  of  Q  then  s'  contributes  three  ^-objects  corresponding  to  the  two  .?' -crossings  and 
the  wall  s' r\int{Q).  The  Q-diagram,  denoted  Vor{Q)  (by  abuse  of  notation),  is  the  Voronoi 
diagram  of  the  ^-objects.  Qearly  the  notion  of  Q-objects  and  2-diagraras  can  be  extended  to 
the  case  where  Q  is  a  union  of  several  quads  in  a  slab.  If  C  is  a  quad  (or  a  union  of  quads), 
we  again  abuse  notation  by  writing  CUaranceQ(p)  for  the  clearance  of  p  with  respect  to  the 
2-objects. 

In  our  algorithm,  a  slab  is  said  to  be  processed  when  the  Q-(^3gram  is  computed  for 
each  active  quad  ^  in  the  slab.  We  represent  a  Voronoi  diagram  as  an  embedded  planar 
graph  ic.,  a  graph  such  that  at  each  vertex  the  cyclic  order  of  the  incident  edges  arc  available 
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and  at  each  face  of  the  embedding,  the  cyclic  order  of  the  bounding  edges  are  available.   For 
each  unbounded  face,  we  introduce  a  fictitious  edge  at  infinity. 

3.   Properties  of  Moving  along  a  Bisector 

Let  s,  s'  be  objects.  In  our  merge  algorithm,  we  need  to  identify  one  of  the  two 
directions  along  the  (j,  j')-bisector  e  as  being  'clockwise'  with  respect  to  s.  This  is  rather 
natural  and  can  be  made  precise  as  follows: 

Let  p  he  a  point  in  e.  For  our  purposes  we  may  assume  that  p  is  not  at  a  transition 
between  two  segments  of  e  corresponding  to  different  governing  equations.  Then  there  is  a 
well-defined  tangent  to  e  at  p.  Let  p,  denote  the  projection  of  p  onto  s.  Let  u  be  a  tangent 
vector  at  p.  First  assume  that  it  is  not  the  case  that  one  of  the  objects  s,  s'  is  a  wall  and  the 
other  an  incident  comer.  Let  v  be  the  vector  from  p.  to  p.  It  is  easy  to  show  that  the  tangent 
at  p  cannot  pass  through  p,,  ie.,  u  x  v  #  0.  Then  we  say  o  represents  a  clockwise  (resp. 
anticlockwise)  direction  about  j  if f  u  X  v  >  0  (resp.  uXv  <  0).  By  symmetry,  it  is  easy  to  see 
that  a  direction  along  e  is  clockwise  about  5  iff  it  is  anticlockwise  about  s' .  Now  suppose  s'  is 
a  wall  and  j  is  an  incident  comer,  and  v  is  a  vector  from  s  directed  away  from  s' .  We  say  u 
represents  a  clockwise  direction  about  s'  iff  it  represents  an  anticlockwise  direction  about  s  iff 
V  X  a  >  0. 

We  now  state  without  proof  some  elementary  properties.  Let  j^,  s\,  s^  be  objects;  e 
and  e'  be  the  (j^,  Sf^)-  and  (s'l^,  j^)-bisectors,  respectively.  Let  p*  =  p*it),  for  r  &  0,  be  a 
parametrized  curve,  thought  of  as  a  moving  point.  Initially,  p*  is  moving  in  e.  Let 
P0  -  PfiiO  (3  =  ^,  ^)  be  the  projection  of  p*  onto  s^. 

(1)  Suppose  that  p*  moves  along  ^  in  the  direction  that  is  clockwise  about  j^.  K  5;j  is  a  wall 
and  pi(  is  in  S/f  then  the  motion  of  p^  (corresponding  to  that  of  p*)  either  is  stationary  at 
an  endpoint  or  is  continuous  and  unidirectional  along  5^.  If  j^j  is  a  comer  (in  which  case 
P/(  is  stationary)  then  the  vector  from  p^  to  p*  is  turning  continuously  anticlockwise 
about  pi^ . 

(2)  Let  ^  be  a  point  in  eCie'  such  that  the  moving  point  p*  meets  q  and  subsequently  moves 
along  e'  in  the  direction  that  clockwise  about  j'^.  We  claim  p*  made  a  left  turn  at  q. 
More  precisely,  if  u  (resp.  v)  is  the  tangent  to  e  (resp.  e')  at  q  in  the  direction  of 
motion  of  p*  then  u  x  v  >  0. 

(3)  Let  I  be  a  separator,  s^  lie  strictly  to  the  left  of  L,  and  j^  lies  on  or  to  the  right  of  L. 
Then  there  is  a  point  on  the  (s^,  s^)-bisectoi  e  beyond  which  moving  clockwise  about  s^ 
has  a  positive  downward  component. 

(4)  The  front  arc  of  p*  is  defined  to  be  the  arc  of  the  clearance  drcle  at  p*  traversed  by 
tracing  the  drcle  anticlockwise  from  p^  to  p,f.  The  back  arc  of  p*  is  defined 
analogously,  as  the  complement  of  the  front  arc  with  respect  to  the  clearance  drcle.  As 
p*  moves  along  e,  this  front  arc  is  sweeping  monotonically  forward  in  a  sense  made 
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precise  in  the  following  lemma: 

Lemma  1.  Let  j^,  s^  (^  X  and  ?  be  a  (j^,  j;;)-edge  in  Vor{X).  If  qQ,  q^  are  the  Voronoi 
vertices  bounding  e  such  that  moving  along  e  from  qQ  to  q^  corresponds  to  clockwise  about  s^ 
then  the  clearance  circle  at  qQ  (w.r.t.  X)  can  only  touch  objects  of  X  on  the  back  arc.  Similarly 
the  clearance  circle  at  <7i  can  only  touch  objects  of  X  in  the  front  arc. 

The  proof  of  this  lemma  uses  a  technique  taken  from  [OSy84a].   It  is  now  convenient  to 
regard  e  as  a  parametrized  curve  p*(r),  for  real  values  t,  with  inaeasing  t  corresponding  to 
the  direction  clockwise  about  s^.    We  will  define  two  dosed  planar  sets  F,,  Z,  of  points  with 
the  following  properties: 
[i]       r,nX,  equals  the  disc  bounded  by  the  clearance  drde  at  p*(l), 

[ii]     the  sets  T,,  Z,  are  continuously  parametrized  (in  the  Hausdorff  metric^  on  sets)  by  t,  and 

[iii]    the  F  (resp.  2)  sets  are  monotonically  growing  (resp.  shrinking).    More  predsely,  for 

f  >  u  we  have  F,  2  F^  and  2„  2  2,.    In  fact,  the  growth  and  shrinkage  have  the 

following  stronger  property:  the  front  (resp.  back)  arc  at  p*(u)  (resp.  p*(t))  lies  in  the 
interior  of  F,  (resp.  2„). 

The  preceding  growth  properties  of  F.  and  2,  immediately  imply  the  lemma.  We  now 
describe  the  F  and  S  sets.  It  is  dearest  to  describe  these  sets  according  to  the  type  of  conic 
that  e  conforms  to.  To  avoid  tedium,  we  will  only  treat  the  case  of  hyperbolas  and  ellipses: 
the  case  where  e  is  parabolic  or  straight  is  essentially  treated  in  [OSY84a].  (We  remark  that 
in  [OSY84a],  one  of  the  cases  also  has  a  third  family  of  sets  A,,  but  this  can  be  avoided.) 


'  The  Hausdorff  mctnc  on  dosed  subsets  of  a  metric  space  K  is  as  follows:  for  any  set  5  C  J'.  «>  0,  let  5, 
be  the  union  of  the  €-balls  about  the  points  of  5.  Then  the  distance  betv.'een  two  closed  iubsets  S,  S'  Q  Y  is  given 
by  d(S,S]  -  infit^  0:  S  QS\andS'  C  5,}. 
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Figure  4.  The  case  of  hyperbolas  with  non-intersecting  drcles 


(a)  Hyperbolas.  If  the  bisector  e  is  hyperbolic  then  both  objects  s^  {^  -  L,  R)  must  be 
circular  walls  (possibly  one  of  the  s^  degenerated  to  a  point).  Let  s^  be  part  of  a  drcle 
Cg  centered  at  q^.  So  the  hyperbola  containing  e  is  (part  of)  the  (Q,  Q)-bisector.  We 
assume  the  centers  are  on  the  jc-axis,  with  <7^  left  of  ^^.  There  are  two  cases,  depending 
on  whether  the  two  circles  intersect.  First  suppose  they  do  not  intersect.  Then  they 
must  lie  external  to  each  other.  G^nsider  the  polygonal  path  P({)  composed  of  the  ray 
from  q^  extending  away  from  q,f,  the  ray  from  q^  extending  in  the  other  direction,  the 
segments  [<7^,  p*(t)]  and  [p*(f),  q^].  Qearly  P(t)  divides  the  plane  into  an  upper  and  a 
lower  part.  Defme  f,  (resp.  I,)  to  be  the  union  of  the  region  above  (resp.  below)  P(t) 
with  the  disc  bounded  by  the  clearance  drcle  at  p*(r).  It  is  easy  to  verify  properties  [i] 
to  [iii]  above.  Now  consider  the  case  where  Q  and  Q  intersect  at  a  pair  of  points  r,  r' . 
Since  j^  and  Sff  are  assumed  not  to  intersea,  e  must  Ue  in  one  of  the  three  sections  of 
the  hyperbola  determined  by  r  and  r' .  If  <r  were  in  the  two  infmite  sections  then  the 
treatment  is  identical  to  the  previous  case.  So  assume  e  lies  between  r  and  r'  inside  the 
'lune'  formed  by  intersecting  the  interior  of  Q  and  of  Q.  The  front  arc  of  p*it)  divides 
the  lune  into  two  parts:  defme  F,  to  be  the  upper  part.  Similarly  defme  1,  to  be  the 
lower  part  in  the  division  of  the  lune  by  the  back  arc  of  p*{t).  The  reader  can  verify 
properties  [i-iii]. 

(b)  Ellipses.  Again,  e  is  elliptic  implies  that  both  objects  s^  are  walls  lying  in  some  circles 
Cp.  However,  we  will  not  assume  that  the  centers  of  the  circles  lie  on  the  jc-axis.  Let  K 
denote  the  strip  bounded  by  the  two  horizontal  lines  through  the  centers  of  Q  and  Q. 
Again  we  consider  two  cases  depending  on  whether  they  intersect.  First  suppose  they  do 
not  intersect.    Then  one  (say  Q)  must  lie  inside  the  other.    Note  that,  because  of  our 
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assumption  that  a  vertical  line  intersects  an  object  at  most  once,  e  must  lie  outside  of  the 
strip  K.  Without  loss  of  generality,  let  e  lie  above  K  and  furthermore  e  lies  in  the 
interior  of  C^  (and  hence  exterior  to  Q).  Let  J'  be  the  interior  of  Q  minus  the  interior 
of  C^.  So  J'  has  a  crescent  shape.  Let  7  be  the  part  of  J'  above  K.  Therefore  e  lies  in 
y.  In  a  natural  way,  the  front  arc  of  p*{t)  divides  J  into  a  left  and  a  right  part.  Define 
r,  to  be  the  part  of  J  to  the  right  part  of  this  division.  Similarly  define  E,  to  be  the  left 
part  of  the  division  using  the  back  arc  of  p*{t).  It  is  again  easy  to  verify  [i-iii]  for  the 
sets  r,  and  2,.  The  other  case,  where  the  two  circles  intersect  should  be  easy  to  defined 
in  analogy  to  the  preceding. 

4.   Merging 

We  now  show  how  to  process  a  slab  5  where  S  is  the  union  of  two  slabs  Si  and  5/j 
separated  by  a  separator  Ly.  By  definition,  this  means  we  compute  the  Q-diagram  of  each 
active  quad  Q  oi  S.  Note  that  Qi  =  Qf^S^'x^  z  union  of  one  or  more  quads  in  5^,  none  of 
which  are  necessarily  active.  The  diagram  of  each  quad  Q'  in  Q^  is  either  already  recursively 
computed  (if  (2 '  is  active)  or  else  the  C -diagram  is  trivial  and  can  be  computed  in  constant 
time.  Thus  with  0{m)  additional  work,  where  m  is  the  number  of  walls  and  comers  in  Qr\Si, 
we  can  assume  all  2' -diagrams  are  available.  Similar  to  the  previous  well-known  algorithms 
for  Voronoi  diagrams  of  curve  segments,  the  Q-diagram  is  obtained  by  'merging'  the  set  of 
C -diagrams  for  bJHQ'  mQi  and  Q^.   We  do  this  in  two  steps: 

(1)  (vertical  merge)  For  (3  =  L,R,  form  the  Qg-diagram  by  'merging'  all  the  Q' -diagrams, 
for  Q'  in  Q^. 

(2)  (horizontal  merge)  Merge  the  Qi-  and  2/?-diagram.    This  is  the  most  important  part  of 
the  algorithm. 

We  can  obtain  the  ^i-diagram  from  all  the  2' -diagrams  (for  Q'  ^  ^z.)  rather  easily  by 
using  the  next  lemma  whose  easy  proof  is  omitted.  Let  Q.,  Qt  ^  ^°  adjacent  quads  of  the 
slab  S^,  and  let  s^^  be  the  long  wall  that  separates  Q  from  Qy  So  the  Lnterseclion  of  s^.^  with 
Si  gives  rise  to  three  objects  5,,  j;,  s^  that  are  simultaneously  Qj-  and  (2;-objects.  Without 
loss  of  generality  let  s^  (resp.  s-^)  be  the  left  (resp.  right)  endpoint  of  Sq. 

Lemma  2.  If  C,(j;)  (i  =  1,  2)  denotes  the  jj-cell  in  the  Q^-di^^am.  then  the  horizontal 
ray  extending  leftward  from  s,  is  contained  in  Ci(j|)nC2(ji). 

Informally,  this  ray  forms  a  natural  boundary  preventing  interaction  of  the  Q-  and  (2;- 
objects.  Using  this  lemma,  it  is  easy  to  justify  the  following  method  for  computing  the  Qo" 
diagram  where  Qq  =  GjUC:'  For  cadi  ^Q-object  s,  the  5-cells  Cq(j)  in  the  Cc^'^agram  is 
obtained  as  follows: 

(1)     Suppose  s  =  J..    The  above  lemma  implies  that   C,{^.)   (i  =  0,  1,  2)   is  unbounded. 
Assume  (?.   lies  above  Q:-    The  boundary  of  the  C„(.j)  is  obtained  partly  from  the 
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boundary  of  Ci(s{)  starting  from  s■^  counterdodcwise  to  its  infinite  edge,  and  partly  from 
the  boundary  of  CiCjj)  starting  from  jj  clockwise  to  its  infinite  edge.   Similarly  if  s  =  Sy 

(2)  If  J  =  S2  then  there  are  two  j-cells:  the  Cq(s)  below  (resp.  above)  j,  is  equal  to  the 
corresponding  cell  in  the  Qj-diagrara  (resp.  Q^-diagtam). 

(3)  If  J  is  not  one  of  j,,  jj.  or  •^3.  then  it  is  a  j2,-object  for  a  unique  /  =  1,  2  and  Cq{s)  is 
equal  to  corresponding  cell  in  the  Qj-diagram. 

By  repeated  application  of  these  observadons,  the  Ci'diagrara  can  be  obtained  in  0{m)  time. 
Similarly  we  can  compute  the  g^-diagram. 

The  main  part  in  constructing  the  ^-diagram  is  the  merging  of  the  Q^-  and  ^^-diagrams. 
This  merging  is  defined  by  a  certain  'merge  curve'  which  generally  consists  of  several 
connected  components.  Indeed,  there  is  a  one-one  coirespondence  between  these  components 
and  the  windows  in  L^fdQ.  Therefore,  we  will  call  the  component  corresponding  to  a  window 
W  the  W -contour,  and  the  procedure  for  computirig  the  W-contour  is  caUed  the  W -merge.  We 
mainly  focus  on  the  mechanism  of  the  procedure  at  present,  leaving  the  correctness  proof  to 
the  next  section. 

Let  W  be  fixed  window  in  QDLf^.  We  wili  assume  that  W  is  finite;  at  the  end  we  will 
handle  the  other  cases.  Let  jq  and  j;  be  the  2 -objects  whose  aossings  at  L^^  determine  the 
upper  and  lower  endpoints  of  W,  respectively.  Let  r,  (i  =  0,  1)  be  the  j,-crossing.  Consider 
the  ray  /?,  emanating  from  r.  downwards  and  normal  to  j..  The  initial  part  of  R,  is  part  of 
the  boundary  of  the  r,-cell  in  the  Q[^-  as  well  as  the  Q/(-(^agram.  Let  rf  (p  =  L,  /?)  be  the 
first  vertex  of  the  Cg-diagram  lying  in  /?,.  Then  define  /?,  to  be  the  closer  of  /t  or  rf  to  r,. 
We  call  P'  the  starter  for  W.  Note  that  starters  is  well-defined  for  a  finite  window  because  the 
ray  /?,  has  a  downward  component  and  must  eventually  get  closer  to  Tq  then  to  r..  Similarly, 
let  /?Q  be  the  ray  emanating  upwards  from  tq  and  normal  to  Sq,  and  define  the  ender  pg  to  be 
the  first  point  along  Rq  that  is  a  vertex  of  the  2^-  or  C^-diagram.  The  W-contour  will  begin  at 
the  top  from  p-  and  terminate  at  p,;. 

Without  loss  of  generality  let  p,  represent  the  intersection  of  the  ray  /?,  with  the 
(s-,  jj-bisector  where  5,  is  a  2^ -object.  It  is  possible  that  Jt  'S  the  jQ-crossing  Tf,.  In  this 
degenerate  case,  the  starter  coincides  with  the  ender  and  the  W-contour  is  defined  to  be 
empty.    Hence  we  may  assume  that  this  is  not  the  case. 
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Ji  is  a  comer  j^  '^  ^  ^^ 

Figure  5.  Elastiating  the  starter 


The  procedure  we  now  describe  is  essentially  that  of  [LD81].    We  construct  the  W- 
contour  as  a  sequence  of  curve  segments 


connecting  a  sequence  of  points 


CTj,  CTj. 


Pl^Pl- 


■<^i- 


•Pi. 


satisfying  the  following: 

(a)  Each  CT;  is  part  of  the  (jf,  jf)-bisector  ?,,  where  sf  {^  =  L,  /?)  is  a  (20-object. 

(b)  The  endpoints  of  ct,  are  p,  and  p;+i  where  the  direction  from  p,  to  p,^,  along  «,  is 
clockwise  about  jf  (and  thus  anticlockwise  about  jf).  (Henceforth,  call  this  the /onvar^ 
direction  along  ?,.  Note  that  this  defmition  of  forward  is  ambiguous  if  {.vf,  if}  =  {t;,,  r. } 
where  r,  is  as  above.  But  this  has  been  excluded.)  Qearly  p,  is  in  the  uuersection  of  the 
bisectors  ^, _ ,  and  e,,  i  >  \. 

Inductively,  suppose  that  p,,  jf,  and  if  have  been  computed.  We  now  extend  this  to 
/  +  1.  Rrst  we  first  show  how  p,^,,  sf^^  are  defined,  returning  to  their  computation  below. 
Let  Cf  be  the  if -cell  such  that  ct,  Q  C^nC^.  Consider  the  component  of  <r,nc3  containing 
p,:  let  (jf  and  rf  be  the  endpoints  of  this  component.  Assume  that  the  forward  direction  along 
f,  is  from  q^  to  rf .  Hence  rf  comes  after  p^  along  <?,.  We  define  p,  +  ,  to  be  the  closer  of  rf- 
and  rf  to  /?,.  Suppose  without  loss  of  generality  that  p^^.,  is  rf.  Then  /7,+,  is  the  intersection 
of  tj  with  a  Voronoi  edge  of  the  Q^-diagram.  Let  if^,  be  the  2^ -object  such  that  this  edge  is 
part  of  the  (sf,  ifv.)- bisector.  Under  the  'general  position'  assumption  that  no  clearance  circle 
intersects  more  than  three  objects,  then  jf-^,  is  uniquely  defined  and  sf^^  will  be  defined  to  be 
jf.     In  this  case,   our  inductive  definition  of  p^.^,  s^^.^    is  complete.     But  in  general,   the 
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clearance  drde  at  /7,+i  (w.r.t.  Q-objects)  may  touch  more  than  three  objects.  By  the  lemma 
in  the  previous  section,  the  objects  can  only  touch  the  clearance  drde  along  the  front  arc 
(recall  that  this  is  the  arc  from  the  contact  point  with  sf  dockwise  to  the  contact  point  with 
sf).  It  is  also  dear  that  if  we  order  the  objects  that  touch  the  front  arc  (starting  with  sf  and 
moving  dockwise  around  the  arc),  we  will  encounter  all  the  ^R-objects  before  any  2i-object. 
Then  jf+j  (resp.  jf+i)  is  defined  to  be  the  last  (resp.  first)  of  these  2;?-objects  (resp.  Q^- 
objects).  This  completes  the  conceptual  description  of  the  W-merge. 

We  still  have  to  initialize  the  merge  and  specify  the  termination  condition.  For  i  =  1 ,  p, 
is  just  the  starter.  Recall  that  above  we  assume  that  the  starter  is  the  intersection  of  the  ray 
J?i  with  the  (jj,  jj^-biseaor,  with  j,  a  G^-object.  In  this  case,  we  let  jf  be  i;  and  jf  be  the 
C/f-wall  jjntnf(2/j).   The  condition  for  termination  is  whenp,+i  is  equal  to  the  ender/^Q. 

We  now  attend  to  the  computational  details  for  obtaining  the  point  p,^, .  It  is  suffident 
to  show  how  to  obtain  the  point  rf .  The  following  is  a  well-known  method  of  obtaining  r^ 
from  q^  where  we  recall  that  q^  is  the  point  where  t,  last  (re-)entered  cell  Cf:  starting  from 
q^,  scan  the  boundary  of  C^  dockwise  (if  3  =  L)  or  antidockwise  (if  P  =  /?)  until  the  first 
time  tfj  intersects  the  boundary.  This  point  is  rf .  We  shall  call  this  the  Shamos-Hoey  scan.  It 
is  not  hard  to  see  that  the  Shamos-Hoey  scan  may  examine  an  edge  more  than  a  constant 
number  of  times. 


17 


W 


Figure  6.  Schematic  illustrating  the  construction  of  the  W-contour 

To  improve  on  this  simple  scheme,  we  must  analyze  the  properties  of  the  points  p,.  Our 
goal  is  to  avoid  rescanning  edges  of  a  cell  Cf  if  the  same  ceU  is  encountered  for  a  succession 
of  values  of  i.  Suppose  that  Cf  =  Cf+j .  The  clearance  circle  at  p^  is  divided  into  the  front  arc 
and  back  arc  by  the  radius  joining  p^  to  the  closest  points  in  the  objects  if  and  jf .  The 
(jf ,  jf)-bisector  e,  continues  beyond  p,  to  meet  the  boundary  of  Cf  at  rf-.  On  the  other  hand 
the  VV-contour  continues  beyond p,  by  following  the  (if+j,  i^'^^J-biseaor  e,^^. 


Definition.  The  truncated  version  Cf  of  Cf  is  defined  to  be  the  part  of  Cf  to  the  right  of 
the  curve  comprising  the  W-contour  from  q^  to  p,  and  the  part  of  e,  from  p,  to  rf.  f^^.^  is 
defined  as  the  first  intersection  of  f|+.  with  the  boundary  of  c'l  scanning  anticlockwise  from 
Pi- 

To  see  that  C,  is  well-defined,  we  must  show  that  the  au^e  in  the  definition  divides  Cf 
into  a  left  and  right  hand  side;  but  this  follows  from  the  fact  that  for  any  pair  of  points  p  =^  p' 
on  this  curve,  the  shortest  open  line  segments  from  p  and  p'  (respectively)  to  s^  are  disjoint. 
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Lemma  3.  In  the  definition  of  the  point  pj+i,  we  could  have  replaced  the  cell  Cf^^  with 
the  truncated  version  Cf,  and  replace  the  point  rf+^  hy  ff+^. 

Proof.  There  are  two  cases  to  consider:  (i)  Suppose  that  after  p^  the  bisector  f,  +  ,  next 
rejoins  e,  at  a  point  q  before  it  reaches  rf+j.  By  defmition,  ff+^  =  q.  But  ^  is  in  the 
(jf,  jf+j)-bisector.  Hence  q  cannot  be  in  the  interior  of  Cf-+y  This  implies  that  the  point  rfv, 
(where  <rf+i  first  meets  the  boundary  of  Cf+J  lies  at  or  before  q.  Thus  p,+i,  if  taken  as  the 
earlier  of  /f+  ^  or  rf+ , ,  would  agree  with  the  original  definition,  (ii)  Suppose  case  (i)  does  not 
hold.  So  either  e,+i  docs  not  intersect  e^  beyond  p,,  or  tf^+j  reaches  rf+i  before  reintersecting 
ej.  The  H'-contour,  in  making  the  transition  from  gj  to  e,+i  at  p,  makes  a  left-ward  turn,  by 
remarks  in  section  3.  Hence  the  point  rf^i  lies  after  rf  (scanning  anticlockwise  from  p,).  In 
this  case  rf+i  is  equal  to  rf+j ,  so  the  lemma  is  again  true,  a 

It  is  important  to  note  that  this  lemma  can  be  iterated  for  each  value  of  i  and  (by 
symmetry)  we  could  similarly  used  C^  in  place  of  Cf.  Thus  we  can  successively  truncate  the 
cells  cf  by  replacing  those  edges  that  have  been  scanned  by  an  appropriate  section  of  e,.  We 
now  summarize  the  scanning  procedure  for  computing  pi+i,  assuming  that  p^,  jf,  and  the 
truncated  cells  Cf  are  available: 

[a]  We  are  at  />,.    Compute  the  values  jf+j  and  jf+j  in  the  manner  outlined  above  (this 
involves  checking  all  Voronoi  edges  that  intersect  p^).   This  gives  us  Cf+^  and  e^+y 

[b]  For  each  ^  -  L,  R,  search  for  the  first  intersection  r,l,  of  tf,  +  .  with  the  boundary  of 
Cf+ ; .   Do  the  search  starting  from  p, ,  scanning  clockwise  if f  (3  =  L . 

[c]  Choose  /7,+,  to  be  the  earlier  of  /f+,  and  r^+^. 

[d]  Update  Cf+ ; .    Replace  the  part  of  its  boundary  that  is  examined  in  [b]  from  p>,  to  rf^ .  by 
the  segment  of  s^+j  between  the  same  pair  of  points. 

This  scan  will  be  called  the  Lee-DrysdaU  scan.  This  completes  our  discussion  for  the 
case  of  a  finite  window  W .    It  remains  to  consider  the  two  cases  when  W  is  infinite: 

(a)  Suppose  W  is  a  half-line.  Without  loss  of  generality,  assume  W  extends  downwards 
from  the  j, -crossing  for  some  Q-object  s^.  The  starter  /?,  for  the  W-contour  can  be  defined 
exactly  as  above,  except  that  it  may  not  always  be  well-defined  (this  happens  if  the  ray  from 
the  J. -crossing  meets  no  Voronoi  vertices  of  Vor{Qi)  and  Vor{Qn).  If  the  starter  is  undefined, 
then  we  say  the  W-contour  is  empty.  If  it  is  defined,  the  merge  process  just  described  can  be 
applied.  Although  there  is  no  ender  we  can  use  the  following  to  give  us  a  criterion  for 
termination:  Let  H  (resp.  //^,  //^)  be  the  convex  hull  of  the  set  of  ^-objects  (resp.  Q^-  ,  Q^- 
objects):  the  hull  is  composed  of  line  segments  and  circular  arcs.  Qearly  there  are  exactly  two 
choices  for  a  pair  pi_,  p^  of  adjacent  vertices  in  H  where  pa  is  a  Qa-<X)mcx.  Choose  the  pair 
that  is  the  lower  of  the  two,  where  'lower'  makes  sense  since  the  segment  \pi,  p^]  intersects 
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Lemma  4.  If  the  W-contour  is  defined  then  the  W-contour  eventually  coincides  with  the 
(Pi,  PRybisector. 

We  omit  the  easy  proof.  The  W-merge  is  done  as  follows:  inductively  assume  the 
availability  of  the  convex  hulls  H^  and  H^.  We  can  compute  H  in  time  linear  in  the  number  of 
comers  in  Q.  Indeed,  this  amounts  to  computing  the  points  p^  and  p^  from  //^  and  H,f.  After 
checking  that  the  starter  is  defined  we  can  do  the  usual  merge,  looking  for  the  (p^,  p^)- 
bisector  as  given  by  the  lemma. 

(b)  Suppose  W  is  the  entire  separator  Lf^.  The  W-contour  is  always  defined  in  this  case. 
It  is  clear  how  the  previous  lemma  leads  to  a  criteria  for  initializing  and  terminating  the  \V- 
merge.  Note  that  this  case  is  just  the  'separable'  situation  that  arises  in  the  original  algorithm 
Hoey  and  Shamos. 

We  have  completed  the  description  of  ihe  W-merge.  It  is  now  fairly  easy  to  construct 
the  Q-diagram  from  the  set  of  W-contours,  Vi^  Q  GHLy.  This  amounts  to  computing  the  cells 
of  each  ^-object:  For  each  aossing  r  at  L^^,  discard  the  r-cells  in  Vor{Q^)  and  Vor{Qf^);  if  r  is 
the  j-crossing  where  j  is  a  Q-waH,  we  can  form  the  two  j<ells  by  merging  the  appropriate 
pair  of  cells  from  Vor(Q[}  and  VoriQ^).  For  the  other  cases,  if  s'  is  a  Q-object,  then  s'  is  a 
Cp-object  for  a  unique  p.  In  this  case,  the  j' -cells  in  Vor(Q)  are  directly  taken  from  Vor(Q^). 
Note  that  if  i'  is  adjacent  to  a  W-contour  then  the  W-merge  would  have  truncated  the  original 
j'-cell  to  the  desired  output. 

5.   Termination 

The  notations  relative  to  a  window  W  from  the  previous  section  is  retained  in  this 
section.  We  now  show  that  the  W-contour  eventually  readies  the  ender  pg  if  W  is  finite.  It  is 
convenient  in  this  section  to  regard  the  W-contour  as  a  parametrized  curve  p*(r),  and  p*  as  a 
moving  point.  Suppose  p*  is  currently  moving  along  cr,  in  the  (jf,  jf)-bisector,  for  some  i.  If 
p*  is  not  at  an  endpoint  (p,  orp,  +  ,)  then  the  projection  of  p*  onto  the  object  s^  (fi  =  L,  R)  h 
well-defined  and  will  be  denoted  p*^  =  p*Jt).  li p*(t)  -  p,^-,  we  (arbitrarily)  define  p* r>  to 
be  the  projection  of  p*  onto  sf  (rather  than  jf+i).  Qearly  p*^{t)  moves  in  a  piecewise 
continuous  manner.  It  is  natural  to  regard  the  path  consisting  of  the  two  segments  \p*i^,  p*] 
and  [p* ,  p*f(]  as  a  'moving  front'.  We  propose  to  study  the  motion  of  the  point  x*  =  x*{t) 
defined  by  the  intersection  of  L^f  with  the  moving  front: 

(j:«}  =  Z.,,n([p«„p*]U[p»,p%])  (•) 

Note  that  x*  is  well-defined  if  p*^  and  p*n  lie  on  opposite  sides  of  the  separator  L^,:  this  will 
be  shown  below.    Recall  that  the  ender  p^  lies  on  the  ray  Rq  originating  from  the  5,,-crossing 

Lemma    5.      Let    the    W-contour    intersect    the    ray    Rq    at    some    point    q.      If 
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Clearance Q{q)  =  d{q,  rj  then  q  =  Pq. 

Proof.    Suppose  the  clearance  circle  C  at  ^  (w.r.t.  Q-objects)  touches  r^.  Let  Cq  be  the 

clearance  circle  at  Pq  (w.r.t.  C-objects).    Note  that  both  circles  are  tangent  to  the  Q-wall 

passing  through  rg.    If  either  circle  properly  contains  the  other,  we  derive  an  immediate 
contradiction.   Hence  Cq  =  C  and  q  =  Pq.   d 


Lemma  6.  The  W-contour  does  not  self-intersect. 

Proof.  If  the  W-contour  self-intersects  then  let  p  be  a  pwint  in  the  interior  of  one  of  the 
fmite  regions  bounded  by  the  W-contour.  Let  p^  be  the  first  point  on  the  W-contour  met  by  a 
horizontal  ray  originating  in  the  'P-ward'  direction.  Let  C^  be  the  clearance  circle  at  p^ 
(w.r.t.  2-objects).  Since  Cg  must  touch  both  a  (7^-  and  a  Cs-object,  it  follows  that  Z,.^ 
intersects  Cg.  This  implies  that  Q  and  C/f  intersect.  Let  L  be  the  vertical  line  through  the 
two  intersection  points  of  Ci^DC^.  If  L\^  is  strictly  left  of  L  then  note  that  Q  must  touch 
some  Q^f-object  at  a  point  in  the  interior  of  Q,  a  contradiction.  By  another  contradiction  in 
the  symmetrica]  case,  we  conclude  L.^  coincides  with  /..  But  in  this  case,  we  easily  see  that  r^ 
and  Tq  are  the  two  points  of  the  intersection  Q  nc^j.  This  implies  that  p^  is  both  the  ender 
and  the  starter,  contradiction,   o 


Lemma  7.  K  p*{t)  is  not  equal  to  p.  or  p^  then  p*i(f)  and  p*r(0  lie  strictly  to  the  left 
and  right  (resp.)  of  Z,,^. 

Proof.    It  follows  from  the  two  previous  lemmas  that  for  each  ^  =  L,  R,  none  of  the 

objects  jf  (i  =  1,  2 )  can  be  equal  to  tq  or  r.  (for  if  sf  =  tq,  say,  then  p*  -  p^,  but  this 

would  imply  that  the  contour  self-intersect).  Hence  each  jf  lies  strictly  to  the  '3-side'  of  Z,,^. 
Now  /7*g(0  lies  in  the  closure  of  sf  and  X  is  proper  implies  p*M)  in  fact  lies  in  if.  Hence  the 
lemma.   □ 

Therefore  the  point  x'(t)  is  well-defmed  for  t  >  0.  We  may  extend  this  definition  by 
defining  x*(0)  to  be  j:*(0*).  Our  goal  is  to  prove  that  x*{t)  is  raonotonically  moving 
downwards  along  Z,,^.  Suppose  the  ender  p^  is  the  intersection  of  the  ray  R^  with  the  (s,  jj- 
bisector  where  Tq  is  the  5r;-crossing  and  5  is  some  Q^-  or  C^j -object.  Let  /j'  be  the  projection  of 
Pq  onto  s  (see  next  figure).  Define  the  point  Xq  to  be  the  intersection  of  L^^  with  the  segment 
|/7(j,  p'].   The  next  lemma  shows  that  the  W-contour  terminates  at  pq  iff  x'{t)  terminates  at  Xq. 


p  =  -r 
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Figiire  7.  The  point  jtq 


Lcmnui  8.  (a)  For  all  r  a  0  corresponding  to  a  point  p*(t)  on  the  W-contour,  x*(r)  is 
above  or  equal  to  xq.   (b)  If  j:*(r)  =  xq  then  p*{t)  =  Pq. 

Proof,  (a)  Let  Cq  be  the  clearance  drcle  at  Pq  (w.r.t.  C-objects).  Let  C*  =  C*(t)  be 
the  clearance  circle  at  p*(f)  (w.r.t.  ^-objects).  If  C*  =  Cq  thenpg  =  p*  and  x*  =  Xq.  Hence 
assume  C*  ^  Cq.  Let  q^  and  <7^  be  the  higher  and  lower  (resp.)  of  the  two  points  of 
intersection  of  L^^  with  C*.  Similarly  let  r^  be  the  point  such  that  Cq  intersects  L^f  at  r^,  and 
Tq.  Note  that  C*(r)  changes  continuously  with  t  and  hence  <7j(f)  moves  continuously  along  L^^. 
It  is  therefore  dear  that  q^(t)  cannot  pass  below  rg.  We  also  claim  that  q^  must  be  higher 
than  T;,:  suppose  otherwise.  Then,  since  C*  cannot  be  properly  contained  in  Cq,  it  is  easy  to 
see  that  both  points  of  the  intersection  of  Cq  and  C*  lie  to  one  side  of  L,^,  say,  to  the  left 
side. 
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Figure  8.  Case  where  ^^  lies  below  r^ 

This  implies  that  the  part  of  C*  lying  right  of  Z,y  is  fully  contained  in  Cq.  But  this  part 
touches  some  2;j-object,  giving  a  contradiction.  Assuming  then  q,,  is  above  r^,,  we  finally 
consider  two  possibilities: 

(i)      Suppose  the  segments  [fq,  r^]  and  [q^,  q^]  are  disjoint.   Then  dearly  x*  is  above  Xq  since 

(ii)  Suppose  [fq,  r^]  and  [^j,  q^]  intersect.  Let  y^^  and  y,f  be  the  two  intersection  points  of  Q 
and  C*.  Qearly  y^  and  _y/j  lie  on  opposite  sides  of  L^^.  Suppose  without  loss  of 
generality,  [p^,  y^]  intersects  L^^  (hence  Pq  must  be  on  the  right  side  of  L^^).  Note  that 
Xq  lies  below  the  intersection  z  of  L,^  with  the  segment  [p^,  y^].  Similarly,  x*  lies  above 
the  intersection  z'  of  L,^  with  the  union  of  segments  (jy^,  p*]U[p*,  >>/?].  But  it  is  easy  to 
verify  that  z  is  below  z'  and  hence  xq  lies  below  x* . 

This  completes  the  proof  of  part  (a).  Part  (h)  foUows  from  the  fact  that  x*(r)  is  strictly 
moving  downwards  with  t:  d  x*{t)  -  xq  and  p*{t)  =^  Pq  then  x*{t  +  t)  is  strictly  below  x^, 
contradicting  (a).   □ 

We  are  now  ready  to  prove  the  main  result  of  this  section: 


Lemma  9.   (Tennination)  The  W-contour  for  a  finite  W  must  terminate  at  the  ender  p^. 

Proof.  Suppose  first  that  p*  (r)  is  not  equal  to  any  p, .  If  j*  (f )  is  the  intersection  of  Z, ,, 
with  \p*i{t),  p*{t)\  then  the  fact  that  p*{t)  is  clockwise  about  si  implies  x*{t)  is  moving 
downwards.  Similarly  if  x*{t)  is  the  intersection  of  L^  with  [/?*(/),  p*/f(0]-  So  suppose  p*{t) 
is  equal  to  some  /?,.  Then  x*  may  suffer  a  discontinuity  at  r.  K  r  is  not  a  discontinuity  for  x* 
then  the  preceding  arguments  show  that  x*  is  moving  downwards  at  r.    So  suppose  there  is  a 
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discontinuity.  Then  x*(r)  =  x*(t)  ¥'  x*(t'*').  But  then  lemma  1  implies  that  x*(r"^)  lies 
strictly  below  x*(t).  Fmally,  lemma  8  implies  that  j:*(f)  is  never  strictly  below  Xq.  So  x*{t) 
eventually  reach  Xq  and  the  same  lemma  asserts  that  this  is  equivalent  to  p*(r)  =  Pq.  o 

6.   Correctness 

For  each  window  W,  the  notations  relative  to  W  from  the  two  previous  sections  will  be 
retained.  From  the  preceding  development,  we  know  that  the  Lee-Drysdale  scan  tenninates 
and  is  correct  (ie.,  equivalent  to  the  Siamos-Hoey  scan)  when  applied  to  a  single  window  W. 
Note  that  the  Lee-Drysdale  scan  is  assumed  whenever  we  refer  to  any  'W-merge'  below.  Now 
we  prove  that  there  is  no  interference  between  the  W- merges  for  different  windows  W.  In 
particular,  we  need  to  be  sure  that  the  modifications  to  the  cells  of  VoriQ^)  made  during  a 
W-merge  will  not  affect  subsequent  merges. 

Lemma  10.  Let  W  and  W  be  two  windows  of  QHL^j.  Assume  W  lies  above  W  and  Tq 
be  the  lower  endpoint  of  W.  Let  Z,  be  the  horizontal  line  through  Tq. 

(!)      The  W-  and  W -contours  are  disjoint.   In  fact  they  lie  on  opposite  sides  of  L. 

(ii)     No  Voronoi  edge  is  examined  in  both  the  W-merge  and  the  W-merge. 

(iii)    The  order  of  doing  the  W-merge  and  the  W  -merge  is  immaterial. 

Proof.  Let  R  be  the  ray  from  rg  extending  horizontally  eastward.  As  in  lemma  2,  we 
know  that  R  is  contained  in  the  Tq-ccU  in  Vor{Q^).  Qearly  (i)  follows  by  sv-rametrical 
considerations  if  we  show  that  the  W-contour  does  not  intersect  R:  this  in  turn  follows  from 
following  fact  (recall  the  termination  proof)  that  for  any  point  p  on  the  W-contour,  if  p  is 
right  of  L^^  then  the  shortest  segment  joining  p  to  a  Ci.-object  intersects  L,^  above  the  point  Xq 
(recall  xq  from  lemma  8).  To  see  (ii),  note  that  any  portion  e  of  an  edge  in  Vor{Qif)  scanned 
in  the  W-merge  has  the  following  property:  it  intersects  or  lies  to  the  left  of  the  W-contour. 
By  'lying  to  the  left'  we  mean  that  for  each  point  p  in  e,  ihc  shortest  segment  joining  p  to  a 
closest  2/f-object  must  intersect  the  W-contour.  Suppose  that  e  is  the  portion  of  an  edge  of 
Vor(Qf^)  that  is  examined  in  both  merges.  K  p  is  any  point  in  e  then  by  the  preceding 
observation,  there  is  a  shortest  segment  from  p  to  a  closest  Q^-objcO.  which  intersects  the  W- 
contour.  Similarly,  there  is  a  shortest  segment  that  intersects  the  W' -contour.  But  these  tsvo 
intersections  lie  on  opposite  sides  of  the  horizontal  line  L.  It  foUows  from  the  triangular 
inequality  that  p  must  in  fact  be  closest  to  the  ^/j-object  tq.  But  the  W-  and  W'  -merges  do  not 
examine  any  edge  on  the  boundary  of  the  rQ-cell,  contradiction.  Part  (iii)  follows  from  (ii) 
and  the  fact  that  Voronoi  edges  are  first  examined  before  being  modified  (discarded  or 
truncated),   a 
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We  finally  show  that  the  g-diagram  constructed  is  indeed  the  Voronoi  diagram  of  the 
set  of  Q-objects.  If  W  is  finite,  define  the  extended  W-contour  to  be  the  union  of  the  W- 
contour  with  the  two  segments  [rg,  p^  and  [rj,  /jj,  where  p^  are  the  starter  and  ender  and  r, 
are  the  two  crossings  bounding  W.  Suppose  W  is  infinite.  K  the  W-contour  is  non-empty,  the 
extended  W-contour  can  be  defmed  analogously.  If  W  is  empty,  then  W  corresponds  to  a 
half-line  determined  by  j-crossing  (for  some  Q-wall  s).  In  this  case,  define  the  extended  W- 
wall  to  be  the  half-line  originating  from  the  j-crossing  in  a  direction  normal  to  s  and  away 
from  the  2-objects.  Defme  C  to  be  the  union  of  all  the  extended  W-contour  where  W  ranges 
over  those  windows  in  L^^  that  intersects  Q,.  We  claim  that  C  divides  the  plane  into  two 
infinite  regions  that  can  be  naturaUy  distinguished  as  the  left  and  right  sides  of  C.  This 
follows  from  the  following  two  facts: 

(a)  C  is  a  doubly  infinite  line. 

(b)  C  does  not  self- intersect.  This  follows  from  the  fact  that  the  individual  W-contours  do 
not  self-intersect  and  distinct  W-contours  are  disjoint. 

L«mma  11.  All  the  2/? -objects  lie  to  the  right  of  C  and  similarly  for  the  2i -objects. 
Only  the  objects  corresponding  to  crossings  at  Ly  lie  in  C  itself. 

Proof.  Let  j  be  a  2;j-object  to  the  left  of  C,  and  assume  that  s  does  not  lie  in  Ly^.  We 
will  derive  a  contradiction.  Let  p  be  the  first  point  on  C  horizontally  to  the  right  of  some 
point  q\Tis.   Qearly  p  is  to  the  right  of  L^^. 

(a)  Suppose  p  is  in  the  W-contour  for  some  W.  The  clearance  circle  C  at  p  (w.r.t.  Q- 
objects)  must  touch  a  point  q  in  some  2z, -object.  Hence  the  radius  of  C  is  at  least  the 
distance  from  p  to  L^^.  Since  the  shortest  line  segment  from  p  to  L,^  is  the  horizontal 
one,  it  follows  that  q  lies  on  this  shortest  segment.  So  <?  is  in  the  interior  of  C, 
contradiction. 

(b)  Suppose  for  some  W,  /?  is  in  the  extended  W-contour  but  not  in  the  W-contour.  Then  p 
is  near  to  a  comer  r  in  L^^.  This  means  the  starter  (similarly  if  it  were  the  ender)  lor  W 
lies  beyond  p  in  the  ray  R,  emanating  from  r.  But  it  is  easy  to  see  that  the  presence  of 
q  implies  that  the  first  vertex  of  Vor(2/j)  lies  between  s  and  p,  contradiction.   □ 

We  now  prove  the  main  result  of  this  section: 

Lemma  12.  (Correctness)  The  Q-diagram  Vor(Q)  is  the  union  of  all  the  W-contours 
together  with  the  portion  of  Vor(Q^)  to  the  left  of  C  and  the  portion  of  Vor{Q,^)  to  the  right 
of  C. 

Proof.  For  each  point  p  in  C,  it  is  easy  to  see  that  p  is  in  some  W-contour  iff  p  is  in 
Vor{Q).   So  assume  p  is  strictly  left  of  C.   The  lemma  then  follows  from  the  foUowing  claim: 
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Clearance Q  (p)  <  Clearance q  (p). 

Suppose  the  claim  is  false  and  for  some  Q^-object  j,  Clearance q  (p)  ^  d(p,  s).    Let  ^  €  5 

such  that  dip,  s)  =  dip,  q).  Then  C  intersects  the  half-open  segment  ip,q\.  Let  r  be  a  point 
in  the  intersection  of  C  with  ip,  q].  Let  Q  be  the  circle  centered  at  p  with  radius  dip,  q)  and 
let  Cj  be  the  clearance  drcle  at  r  (w.r.t.  (2-objects).  The  interior  of  Q  does  not  intersect  any 
C^-object.  Qearly  Cj  is  contained  in  Cq  (otherwise  q  would  be  in  the  interior  of  C^). 
Consider  two  cases: 

(a)  q  is  not  a  comer  in  Ly.   Then  s  is  not  a  Q^-object.   Since  Cj  must  touch  some  C^-object 
s' ,  it  follows  that  the  interior  of  Cq  intersects  s' ,  contradiction. 

(b)  <7  is  a  comer  in  L^.   Let  q  be  the  j' -crossing  where  s'  is  a  Q-wall.   Thus/?  is  on  the  line 
through  q  and  normal  to  5' .   Thus  p  lies  in  C,  contradiction,   o 


7.   Putting  it  together 

The  main  procediu-e  consists  of  two  preprocessing  steps  followed  by  a  call  to  a 
recursive  procedure: 

Main  Procedure 

Input:  a  firoper  set  X  of  objects. 

(1)  Presort:  Sort  the  set  of  comers  according  to  their  vertical  projection  onto  the  x-axis. 

(2)  Pre-scan:  Do  a  scan-line  sweep  of  the  line  segments  to  determine  for  each  comer  s  the 
walls  that  are  immediately  (vertically)  above  and  below  5.  This  is  essentially  the 
algorithm  of  Hoey-Shamos  for  detecting  line  segment  intersections.  Use  aboveis)  and 
below{s)  to  represent  these  walls.  If  s  has  no  walls  above  or  below  it  then  this  is  given  a 
special  indicator.  This  step  uses  the  information  gathered  in  the  pre-sorting.  During  the 
pre-scan  we  also  introduce  the  set  of  separators.  Let  left{s)  and  nght{s)  indicate  the  two 
separators  adjacent  to  s.  (Note:  to  handle  the  case  of  more  than  one  comer  in  any 
vertical  line,  only  trivial  modifications  are  necessary  in  our  entire  preceding 
development.) 

(3)  Recursion:  Call  a  recursive  procedure  to  process  the  slab  S  bounded  by  the  leftmost 
and  rightmost  separators.  Note  that  the  entire  slab  5  constitutes  an  active  quad,  so 
the  diagram  of  this  quad  is  VoriX). 

End  Main  Procedure. 

We  now  present  the  recursive  procedure  for  processing  an  arbitrary  slab  S  (represented 
by  a  pair  of  separators).  We  assume  that  X,  as  well  as  the  other  information  of  the  pre- 
scanning,  is  available  to  the  procedure  via  global  variables.    The  recursive  procedure  returns 
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(i)  a  list  of  the  active  quads  of  S  and  a  list  of  windows  for  each  quad,  (ii)  a  list  of  the  Voronoi 
diagrams  of  these  quads,  and  (iii)  the  convex  hull  of  the  Q-objects  for  the  topmost  and 
bottommost  active  quads  Q  in  the  slab. 

Recursive  Procedure 

Input:  (m,  S)  where  5  is  a  slab  and  m  the  number  of  comers  in  S. 

(1)  Basis:  If  m  =  1,  then  use  the  pre-scanning  information  to  determine  the  unique  comer  p 
in  S  and  also  the  quad  Q  containing  p.  Compute  the  diagram  of  Q  and  retxim.  (Note 
that  the  division  of  walls  into  subobjects  at  their  aossings  is  done  at  this  point  in  the 
algorithm.) 

(2)  Divide:  Divide  the  slab  S  into  two  slabs  5^  and  S/f  with  \m/2\  and  [m/2j  vertices 
respectively.  (This  is  easy  to  do  assuming  that  the  separators  for  S  are  just  indices  into 
an  array).   Recurse  on  S^  and  S^.   Let  L,^  denote  the  separator  between  5^^  and  S^. 

(3)  Conquer: 

(3.1)  Vertical  Merge.  It  is  actually  quite  interesting  to  determine  the  list  of  active  quads  Q  in 
S.  In  order  to  preserve  the  continuity  of  the  main  description,  we  defer  this  to  the  end 
of  this  section.  For  now  we  just  note  that  for  each  active  Q  Q  S  with  k  objects,  we  can 
gather  together  the  set  of  active  and  inactive  quads  oi  S^(^  -  L,  R)  whose  union  forms 
Q.  Furthermore,  if  5  has  <:  active  objects  then  this  takes  time  0(k)  over  all  the  active 
Q.  We  then  merge  the  diagrams  of  the  quads  in  Q^  -  QHS^  to  obtain  the  diagram  of 
Qn,  using  the  method  described  in  section  4. 

(3.2)  Horizontal  Merge.  We  assume  that  a  list  of  the  active  quads  Q  of  5  and  also  a  list  of  the 
windows  of  each  Q  is  available  from  the  previous  step.  We  apply  the  W-merge  to  each 
window  W  C  L,\{<^Q,  for  each  Q.  Fmally,  we  construct  the  Q-dlagram  from  these. 
Note  that  in  the  process  of  doing  the  W-merge  for  the  topmose  and  bottommost  quads, 
we  have  also  computed  their  convex  hulls. 

End  Recursive  Procedure. 

We  are  now  done  except  for  the  details  on  determining  the  active  quads  of  5  from  the 
corresponding  lists  for  5^  and  5;^.  Let  us  introduce  some  terminology  relative  to  any  slab  S 
and  separator  L:  an  S-interval  /  in  Z.  is  one  of  the  form  /  -  Lr\Q  for  some  quad  Q  m  S.  An 
S-active  interval  is  similarly  defined  except  that  the  quad  under  consideration  must  be  active. 
Two  intervals  in  L  overlap  if  their  interiors  intersect.  Given  a  set  /  of  intervals  in  L,  the  /- 
equivalence  relation  on  /  is  defmed  as  the  reflexive,  symmetric,  and  transitive  closure  of  the 
overlap  relation  on  /.  Let  the  slab  S  be  divided  into  slabs  Si  and  S^  by  the  separator  L^f.  Let 
Jg  be  the  set  of  5g-active  intervals  in  L,^. 

Lemnu  13.  There  is  a  bijective  correspondence  between  the  set  of  5-active  intervals  in 
L\^    and    the    set    of    (7^ U //j)-equivalence   classes   such   that    if   7   is   an  5-active   interval 
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corresponding  to  an  eqxiivalence  class  then  /  is  equal  to  the  union  of  the  intervals  in  that 
equivalence  class. 

Proof.   Suppose  /  is  an  5-active  interval.   The  set 

c,  =  {y  e  j^uj/f  :yc/} 

of  Si'  and  5;j-active  intervals  inside  /  must  be  non-empty.  Note  that  an  interval  not  in  G; 
cannot  be  ( J^^  U /;j)-equivalent  to  any  interval  in  C;.  Hence  G;  is  a  union  of  equivalent  classes. 
Qaim:  there  is  only  one  eqviivalent  class  in  G;.    Let  £;  be  the  set  of  endpoints  of  intervals  in 

Gj.    Sort  Ej  as  jti,  jti x^,  {k  a  2)  according  to  their  height.    If  it  =  2  then  the  claim  is 

immediate,  so  let  Jk  >  2.  Note  that  if  x^  is  not  one  of  the  endpoints  of  /  (ie.,  1  <  j  <  it)  then 
Xj  cannot  be  the  endpoint  of  both  an  5^-  and  an  5^-interval.  Suppose  Xt  is  the  endpoint  of 
some  5^-interval.  Then  we  see  that  or,  is  in  the  interior  of  some  5;{-active  interval  /•:  this  is 
because  jc,  represents  an  j-crossing  for  some  C-object  s  and  one  of  the  endpoints  of  5  is  in  L/j, 
making  the  S/j-interval  containing  jc,  active.  So  let  x,  (j  >  2)  be  the  lower  endpoint  of  J  (the 
upper  endpoint  is  x,).  A  similar  argument  shows  that  x^  is  in  the  interior  of  some  5^-active 
interval  /,.  Furthermore,  /j  and  I2  overlaps.  Continuing  this  way,  we  eventtially  get  to  an 
interval  whose  lower  endpoint  is  x^.  □ 

Using  this  lemma,  we  can  compute  the  list  of  the  5-active  intervals  at  each  of  the  two 
separators  bounding  S:  assume  that  inductively  we  have  a  list  of  the  5g-active  intervals  at  each 
of  the  separators  bounding  Sq.  We  can  now  'merge'  the  list  of  5^  and  S^  at  Z,^^  in  the 
straightforward  way.  When  done,  we  can  easily  convert  this  list  into  the  corresponding  5- 
active  list  of  intervals  at  the  two  boundary  separators  of  5. 

8.   Complexity 

We  first  analyse  the  complexity  of  the  W-raerge. 

Lemma  14.    The  W-merge  takes  time  0{k  +  l)  where  k  is  the  number  of  segments 

CTi.cT, cTj   comprising  the   W-contour  and  /  is  the  number   of  examined  edges  from 

Vor(ej  and  VoriQ^). 

Proof.  The  work  to  construct  cr,  is  proportional  to  the  number  of  edges  of  the  cells  Cf 
(3  =  L,  R)  examined  in  the  scan.    Refer  to  section  4,  steps  [a-d]  of  the  scanning  procedure: 

Step  [a]: 

Here  we  determine  all  those  Voronoi  edges  that  terminates  at  p, .  It  is  not  hard  to  show 
that  none  of  these  edges,  except  for  the  edges  that  bound  the  cells  Cf^  • ,  C^. . ,  will  ever 
be  examined  again.  Hence  the  work  done  in  this  step,  summed  over  the  entire  W- 
merge,  is  0(k)  for  those  edges  that  wiU  not  be  examined  again  and  0(1)  for  the  others. 
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Step  [b]: 

Among  the  edges  of  Cf+i  that  are  scanned  in  step  [b],  all  except  possibly  the  first  and 
second  edges  are  scanned  for  the  first  time.  We  can  charge  the  cost  of  scanning  the  first 
two  edges  to  a^+i,  and  the  rest  are  charged  to  the  edges  themselves.  The  overall 
charges  for  these  are  0(1)  and  0{k),  respectively. 

Steps  [c,d]: 

These  can  be  made  0(1)  for  each  segment  a^,  and  hence  0(1)  when  summed  over  the 
entire  W-raerge.   d 


Lemma  15.  The  construction  of  a  Q-diagram  is  0(m)  if  there  are  m  2-objects. 

Proof.  Recall  that  before  the  'horizontal'  merging  is  applied  to  Vor(Q^)  and  Vor(Qf(), 
the  'vertical'  merging  is  done.  The  time  bound  is  0(m).  To  bound  the  time  for  horizontal 
merging,  recall  that  In  the  last  section  we  have  essentially  showed  that  the  separate  W- merges 
do  not  interact.  From  the  previous  lemma  the  work  done  in  the  W-merge  is  0(kv/+lw)  where 
ky^  is  the  number  of  segments  comprising  the  W-contour  and  /«,  is  the  number  of  examined 

edges.    The  sum  2*w  (over  all  windows  W  of  C)  is  order  of  the  number  of  segments  in  the 

w 

curve  C,  and  the  sum  2'iv  is  order  of  the  number  of  Voronoi  edges  in  Vor(Qi)  and  Vor(2s)- 

w 

Since  Vor(Q)  has  no  more  edges  than  the  combined  number  of  edges  in  Vor(Q^),  Vor(Qff)  and 

C,  it  is  dear  that  the  sura  '^(kw  +  lw)  o\eT  all  windows  W  is  0(m)  where  m  is  the  number  of 

w 

C-objects.   □ 


Theorem  16.  The  Voronoi  diagram  of  a  set  X  o(  n  pairwise  disjoint  objects  cab  be 
computed  in  0(n\ogn). 

Proof.  As  in  the  introduction,  we  can  divide  the  work  done  in  the  various  recursive  calls 
into  stages,  where  each  stage  corresponds  to  processing  a  collection  of  slabs  5.,  5;,  .  .  .  ,5„, 
where  these  slabs  are  disjoint  and  the  union  of  the  slabs  contains  all  the  objects  in  X.  Since 
each  object  occurs  in  at  most  two  active  quads  in  these  slabs,  and  the  processing  of  each  quad 
is  linear  in  the  number  of  objects  in  it,  we  conclude  that  a  stage  takes  linear  time.  But  there 
are  0(log/i)  stages,   a 


9.    Diagram  for  [ntersecUng  clrdes 

[to  be  filled  in] 
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10.   Coodusion 

This  paper  solves  the  open  problem  of  an  O(nlogn)  algorithm  for  the  problem  of 
computing  the  Voronoi  diagram  of  a  set  of  points,  line  segments  and  drcuJar  arcs.  The 
algorithm  is  simple  enough  to  have  an  impact  on  practical  applications,  and  in  particular,  in 
motion-planning  for  robots.  The  technique  clearly  works  for  more  general  algebraic  curves, 
provided  we  take  care  to  break  up  each  curve  into  a  number  of  suitably  small  sections.  It  also 
seems  that  the  technique  can  be  generalized  to  handle  the  Voronoi  diagram  of  a  set  of 
polyhedral  objects.  This  is  a  subject  of  further  research. 
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